subroutine r2abij_t2efij_iefab_gin(r2, t2, &
              iabcd, icore, fact)
!
use mod_iop
use mod_size
use mod_orbit
use mod_ioff 
implicit none
!
real*8, intent(in),device :: iabcd(*)
real*8, intent(in),device :: r2(*)
!
real*8, intent(inout) :: t2(*), icore(*)
real*8, intent(in) :: fact
!
integer :: i, j, a, b, nsize, i0, i1
real*8, external :: dnrm2
!this is to do r2 = r2 + fact*t2(efij)*iefab
if(prec==8) then 
   call r2abij_t2efij_iefab_dgin(r2, t2, &
              iabcd, icore, fact)
elseif(prec==4) then 
   call r2abij_t2efij_iefab_sgin(r2, t2, &
              iabcd, icore, fact)
endif
return 
end

